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Abstract. We consider the Cauchy problem for an energy supercritical nonlinear wave equation 
that arises in (1 + 5) -dimensional Yang-Mills theory. A certain self-similar solution Wo of this 
model is conjectured to act as an attractor for generic large data evolutions. Assuming mode 
stability of Wo, we prove a weak version of this conjecture, namely that the self-similar solution Wo 
is (nonlinearly) stable. Phrased differently, we prove that mode stability of Wo implies its nonlinear 
stability. The fact that this statement is not vacuous follows from careful numerical work by Bizoh 
and Chmaj that verifies the mode stability of Wo beyond reasonable doubt. 

1. Introduction 

This paper is concerned with the study of the Cauchy problem for the semihnear wave equation 
Mt, r) - Vrr(i, r) - ^^Mt, r) + r)) = (1.1) 

for (i = 5 where F{ij)) = + 1)(V' + 2) and r = x G M'^. Eq. (II. ip arises in S'0(o?)-equivariant 
Yang-Mills theory, see [15) . [5j for a derivation. Historically, the introduction of nonabelian gauge 
theory by Yang and Mills was fundamental for the development of the standard model of particle 
physics, see [39]. Apart from that, the Yang-Mills model as a classical field theory attracted a 
lot of interest too, cf. [1]. Furthermore, Yang-Mills equations have been proposed as toy models 
for Einstein's equations of general relativity, see e.g. ^16j, [20]. Especially in this context the 
development of singularities in finite time is of interest. 

Eq. (|l.ip is invariant under scaling and it admits a conserved energy 



d-3 



dr 



which scales like 

E{i;\Xt, •), V't^(At, •)) = A'^-^S(^(t, ■),Mt, •)) 
where ip'^{t,r) := ^/'(t/A, r/A), A > 0. This shows that Eq. (jl.ip is energy subcritical in the physical 
dimension d = 3 whereas it is critical for d = 4 and supercritical for d > 5. According to the 
usual blowup heuristics, where energy conservation prevents the solution from shrinking in the 
energy subcritical case, one expects global existence for d = 3. Indeed, Eardley and Moncrief 
[16| . |17] considered this problem (without symmetry assumptions) and proved global existence for 
data in suitable Sobolev spaces without restriction on size. This classical result was strengthened 
by Klainerman and Machedon |27| who lowered the required degree of regularity of the data and 
provided a different approach to the problem, see also [22j, [40]. 

Yang-Mills fields in dimension d = 4 attracted a lot of interest in the recent past due to new 
developments in the study of energy critical wave equations. Local well-posedness under minimal 
regularity assumptions was considered by Klainerman, Tataru |28] and global existence for small 
data was proved by Sterbenz [38] . In the critical dimension d = 4, Eq. (jl.ip admits a static, finite 
energy solution which is known as the instanton. This indicates the existence of more complex 
dynamics than in the case d = 3. In particular, the energy of the instanton represents a threshold 
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for global existence. Indeed, Cote, Kenig and Merle [10] proved global existence and scattering 
(either to zero or to a rescaling of the instanton) for data with energy below (or equal to) the 
energy of the instanton. On the other hand, it is known that the Yang-Mills system in the critical 
dimension d = A can develop singularities in finite time. This has been conjectured by Bizoh [5] 
and demonstrated numerically in [5], [31]. Furthermore, the blowup rate was derived by Bizoh, 
Ovchinnikov and Sigal |3]. The existence of blowup solutions was proved rigorously by Krieger, 
Schlag, Tataru [2^ as well as Raphael and Rodnianski |35] who also obtained the stable blowup 
rate. 

In the supercritical dimension d = 5, which shall concern us here, much less is known. In general, 
the study of energy supercritical wave equations is still only at the beginning, see, however, e.g. [21], 
[23] . [26] . [25] . [8], [7], [l2] for recent progress. It is clear that energy supercritical problems will have 
to play a prominent role in the future development of the field, not only because of their relevance 
in physics which can hardly be overstressed. In particular, much work remains to be done in order 
to improve our understanding of large solutions which at the moment is mostly confined to the 
construction of self-similar solutions by solving a corresponding elliptic ODE problem, a procedure 
which is insensitive to the criticality class of the problem, see e.g. [36], [9], [4], [6]. In the case of 
the Yang-Mills field in d = 5 one has global existence [37j for small data, see also [30] for d > 6, 
whereas for large data finite time blowup is possible [9]. The failure of global existence has been 
demonstrated by constructing self-similar solutions to Eq. (|l.ip . In fact, Bizoh [5] showed that 
there exists a countable family of self-similar solutions. Furthermore, he was even able to find an 
explicit expression for the "ground state" of this family which we denote by ijj'^ and it reads 

i;^ {t,r) = Wo{ 

where T > is a constant (the blowup time) and 

The relevance of such an explicit solution for understanding the dynamics of the equation depends 
on its stability. In other words, the important question is: does the blowup described by tp'^ occur 
for a "large" set of initial data? Numerical simulations [2] indicate that this is indeed the case. 
Moreover, it appears that the blowup via ip'^ is "generic", i.e., sufficiently large, "randomly" chosen 
initial data lead to an evolution which asymptotically (as t — )• T— ) converges to xp"^ . Consequently, 
in the present paper we study the stability of tp'^ and obtain the following result, see Theorem 11.31 
below for the precise statement. 

Theorem 1.1 (main result, qualitative version). Suppose iJj^ is mode stable. Then there exists an 
open set (in a topology strictly stronger than the energy) of initial data for Eq. (11. ip such that the 
corresponding time evolution approaches ij}^ and blows up. 

We remark that the technical assumption we have to make, the mode stability of '0"^, is equivalent 
to a certain spectral property of a (nonself-adjoint) second order ordinary differential operator, see 
below. Unfortunately, we are not able to prove this spectral property rigorously but it has been 
verified numerically beyond reasonable doubt [2j. In the following we will comment on this issue in 
more detail. 

1.1. The mode stability problem. A first step in the stability analysis of the self-similar solution 
ip'^ is to look for unstable mode solutions. To this end we insert the ansatz ip = ip'^ + 93 into the 
Yang-Mills equation (jl.ip with d = 5 and linearize by dropping all nonlinear terms in 93. This 
yields the evolution equation 

iftt - ^tt - ^^r + ;|F'(M^o(t^) - = (1.2) 



for the perturbation Lp. In order to obtain a time independent potential we introduce p = 
as a new variable and restrict ourselves to the backward lightcone of the blowup point (T, 0) by 
assuming /? G [0, 1]. A particularly convenient choice for a new time variable is r = — log(T — t). 
The coordinates (r, p) are sometimes referred to as "similarity variables" and they are frequently 
used when dealing with self-similar solutions for nonlinear wave equations, see e.g. |12] . |13) . 
[32] , [33l, f33]. Note that the blowup takes place as r ^ oo and thus, we are effectively dealing with 
an asymptotic stability problem. By setting (p{t,r) = (/)(— log(T — t),-!^—^^, Eq. (|1.2p transforms 
into 



+ (l)r + 2p(l)rp - (1 - p^) 



+ j,F\Wo{p) 



0. 



(1.3) 



A solution of Eq. ()1.3p of the form (px{T, p) = e^'^u\{p) for A G C and a nonzero function 
u\ G C°°[0, 1] is called a mode solution. It will become clear below why we can restrict ourselves 
to smooth u\. Furthermore, we say that A is an eigenvalue (of ip'^) if there exists a corresponding 
mode solution (l)x- For obvious reasons a mode solution (or an eigenvalue A) is called stable if 
ReA < and unstable otherwise. At this point it is worth emphasizing that a priori the nonexistence 
of unstable mode solutions is neither necessary nor sufficient for the (nonlinear) stability of ijj'^ . 
However, the nonexistence of unstable mode solutions is obviously necessary for the linear stability 
of and, as we will prove in this paper, it is even sufficient for the nonlinear stability of ip'^ 
(see Theorem 11.31 for the precise statement). As a consequence, it is crucial to understand mode 
solutions and, by inserting the ansatz (j)xiT,p) = e^'^u\{p) into Eq. ()1.3p . this problem reduces to 
the ODE 



(1 



+ 2Xpu\ + A(A + l)ux + jzF'{Wo{p) 



l)ux = 0. 



(1.4) 



As a matter of fact, there exists an unstable mode solution for A = 1 given by ui{p) = pWq{p). 
However, it turns out that this is a symmetry mode, i.e., it stems from the time translation symme- 
try of Eq. (jl.ip and does not count as a "real" instability (this will become much clearer in Section 
3] below, see in particular Lemma [^77|) . Consequently, we define 



Definition 1.2. The solution ^/^^ is said to be mode stable iff ui{p) 
of Eq. (JLII) in C~[0, 1] with ReA > 0. 



pWq{p) is the only solution 



Unfortunately, it appears to be extremely difficult to exclude unstable mode solutions (one has to 
bear in mind that the problem is nonself- adjoint; so in principle there could be unstable eigenvalues 
with nonzero imaginary parts). If ReA > 1, the problem is fairly easy since one can resort to Sturm- 
Liouville oscillation theory [5] and it is well-known that there do not exist unstable eigenvalues A 
with Re A > 1 apart from the aforementioned symmetry mode. However, the domain < ReA < 1 
seems to be very challenging. The fact that exactly the same problem occurs in the study of energy 
supercritical wave maps |14j . [12j underlines the importance of having a general approach to that 
kind of nonself-adjoint spectral problems. For the moment, however, this question remains open. 
On the other hand, there are very reliable numerical techniques to study boundary value problems 
of the type p.4p . As a consequence, the mode stability of ij^'^ has been established numerically 
beyond reasonable doubt [2]. In addition, we provide a new result (see Lemma [3.101 below) which 
excludes unstable eigenvalues that are far away from the real axis. This puts the available numerics 
on an even stronger footing. 
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1.2. The main result. With these technical preparations at hand we can formulate our main 
result. To begin with, we define a norm 

t-R 2 

+ / \rg'\r) + bg'{r) + ^^g{r)\ dr 
Jo 

on i{R) := {{f,g) G C^[0,R] x C^[0,R] : /(O) = /'(O) = g{0) = 0}. It is easily seen that || • Wg^^R) is 

indeed a norm on this space and we denote by £{R) the completion of £{R) with respect to || • \\£(r)- 
Our main result is the following. 

Theorem 1.3 (stable self-similar blowup). Assume to he mode stable. Let e > and suppose 
we are given initial data {f,g) € ^(|) such that 

||(/,5)-(V'(0,-),^hO,-)ll£(|) 
is sufficiently small. Then there exists a unique solution tp of Eq. (jl.ip with d = 5 satisfying 

V,(0,r) =/(r), ^t{0,r) = g{r), r G [0, |] 
and a blowup time T G (^, |) such that 

(T - tf^mt, ■),Mt, •)) - (^^(t, ■),i^I{t, -meiT^t) < Ce{T - tt'^^-^ (1.5) 

for all t G [0,T) with ujq := max{ — |,;Uo} where /xq is the real part of the first stable eigenvalue of 
ijp- and Ce > is a constant which depends on e. 

Remarks: 

• As usual, by a "solution" we mean a function which satisfies the equation in the sense of 
Duhamel. 

• By a simple scaling argument one immediately sees that 

ii(v^(t,-),^ni,-))ikT-*)^(r-ri 

for all t G [0,r). This explains the normalization in p.5|) and shows that "0 converges to 
the self-similar solution ifF in the backward lightcone of the blowup point. 

• The rate of convergence in (jl.Sp is dictated by the first stable eigenvalue which complies 
with heuristic expectations and numerics The e~loss in our estimate is purely technical. 
We also remark that the numerically obtained value for is approximately —0.59 [2]. 

• The f-norm is very natural since it is derived from a conserved quantity of a suitable "free" 
equation which is associated to Eq. (jl.ip . see below. Furthermore, the boundary conditions 
assumed for the initial data (/, g) are natural too, since any sufficiently regular solution ifi 
of Eq. (jl.ip must satisfy ^(t, 0) = iptit^^) = ipr{t,0) = for all t (provided -0 belongs to 
the same topological sector as 'ip'^ which we obviously assume). 

• It is clear that the result cannot hold in the energy topology. This is due to the fact that 
the corresponding local energy in the backward lightcone of the blowup point decays like 
T — t as t T— and thus, self-similar blowup is invisible in the energy topology. This is, 
of course, nothing but a manifestation of energy supercriticality. 

1.3. An outline of the proof. The proof consists of a perturbative construction around ip'^ which 
proceeds in several steps. 

(1) First, we identify a suitable Hilbert space where the corresponding inner product is 

• derived from a conserved quantity of a suitable "free" equation where the latter is 
(roughly speaking) obtained from Eq. (II. ip by dropping the nonlinear term, 
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• strong enough to detect self-similar blowup, 

• strong enough to control the nonlinearity. 

(2) Next, we introduce similarity coordinates, linearize Eq. (jl.ip around V'^ and construct a 
semigroup that governs the linearized evolution. The application of semigroup theory to 
this problem is natural since 

• the involved differential operator is highly nonself-adjoint due to the introduction of 
nonorthogonal coordinates, 

• the evolution problem is restricted to the backward lightcone of the blowup point and 
is thus only well-posed in forward time which is reminiscent of parabolic equations. 

(3) Then we perform a detailed spectral analysis of the semigroup generator and construct a 
Riesz projection of rank 1 which removes the unstable symmetry mode that results from 
the time translation invariance of Eq. (II. ip . As a consequence, we obtain exponential decay 
of the linear evolution on the codimension 1 stable subspace. 

(4) Next, we prove a Lipschitz property of the nonlinearity which allows us to run a fixed point 
argument in order to obtain a solution to the nonlinear problem. That this is possible is 
not surprising since our norm controls sufficiently many derivatives to obtain a Moser-type 
estimate. Furthermore, the linear evolution decays exponentially and this kind of decay 
is reproduced by the Duhamel formula. However, in order to suppress the instability of 
the linear evolution we have to modify the data. This is similar to the Lyapunov-Perron 
method in dynamical systems theory. 

(5) Finally, we show that the aforementioned modification of the data is equivalent to changing 
the blowup time. Thus, by choosing the appropriate blowup time we obtain a solution of 
the original equation (jl.ip with the properties stated in Theorem 11.31 

1.4. Notation. As usual, we write a < 6 if a < c6 for some c > 0. Similarly, we use > and a ~ 6 
means a ^ b and b ^ a. The big-0 symbol has its standard meaning from asymptotic analysis. 
In order to improve readability we use boldface letters for vectors and number the individual 
components by subscripts, e.g. u = (ui,U2). The symbol denotes the Frechet derivative and 
Dj^ is used for the j-th partial Frechet derivative. For a closed linear operator A we write cr{A), 
crp{A), Ra{^) for the spectrum, point spectrum and resolvent of A, respectively. 

2. Transformation to a first-order system and similarity coordinates 

2.1. Formulation of the Cauchy problem. As explained in the introduction, we intend to study 
the Cauchy problem Q 

/ ^tt - Ar - f V'r + ^F{lP) =0 , . 

\V'(o,.) = /, Mo,-)=g ^ ^ ^ 

for a function ^ : Ct — > K where 

Ct ■■= {it,r) : t e [0,T),r e [0,T - t]}, T>0 

and f,g : [0, T] —?■ M are prescribed, sufficiently regular functions (the initial data). Furthermore, 
the nonlinearity F is given by F{'ip) = 'i/'(V' + l)(V' + 2)- Note also that the requirement of regularity 
at the center demands tpr{t,0) = and ip{t,0) G {0,-1,-2} for all t. We are interested in the 
stability of the blowup solution 

where 



In order to avoid notational clutter we usually omit the arguments and write ip instead of ip{t,r). 
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is the Bizon solution. Since V'"^(i>0) = we restrict ourselves to solutions of Eq. (j2.6p that satisfy 
tp{t, 0) = for all t. In the following we perform some formal manipulations to transform (j2.6p into 
a convenient form suitable for further analysis. 

We intend to study small perturbations of ip'^ and thus, it is reasonable to reformulate (j2.6p 
relative to tp'^ , i.e., we insert the ansatz tp = tp'^ + ip into Eq. (|2.6p and obtain the Cauchy problem 

r ^,,-^^^-l^^ + ^,p + ^[F'{^pT)-2]ip + j^NT{ip) = , . 

\ (^(0,.) = /-^^(0,-), ¥^t(0,-)=5-^nO,-) ^ ^ 

for the perturbation 99 : Ct — IK. Here, 

Nt{v) = 3[F(V^^ + ^)- - F'ii,'^)^] (2.8) 

is the nonlinear remainder. Observe further that 

F'(V'^(t,0))=F'(0) = 2 

and thus, by subtracting the constant 2 we have regularized the "potential term" in such a way 
that ^ l^2^^^~'^ remains bounded as r — )• 0+. Finally, the perturbation ip inherits the boundary 
conditions V3(t,0) = (pr{t,Q) = for all t. So far nothing has happened and Eq. (j2.7p is equivalent 
to Eq. (j2.6p \i 'p = 'p'^ + 99. In order to fix terminology we call 

the free equation, 

the linear or linearized equation and, finally, the full problem Eq. (12. 7p is referred to as the nonlinear 
equation. Note carefully that we have assigned all singular terms to the free equation. This is 
necessary since our overall strategy is to treat the nonlinear equation as a perturbation of the 
linearized equation which, in turn, is viewed as a perturbation of the free equation. Therefore, the 
topology is dictated by the free equation. 

2.2. Higher energy norm. As already outlined in the introduction, the energy topology is too 
weak to study self-similar blowup. Consequently, we have to find a stronger norm and it is ad- 
vantageous if this norm is naturally associated to the free equation. Furthermore, we intend to 
control the nonlinearity by a Moser-type estimate and therefore, we expect to need at least |+ 
derivatives (recall that the Yang-Mills problem is in 1 + 5 dimensions). For simplicity, however, we 
avoid fractional Sobolev spaces and aim for a norm that controls 3 derivatives. The key observation 
in this respect is that, if we set 

^{t,r) :=ia,[i5,(rVi,0)], 

we obtain the identity 

'Ptt - ^r-r = {^dr [r^ {(ftt " ^rr " ^^r + ^V^)] } • 

Thus, if 99 satisfies the free equation then p is a solution to the one-dimensional wave equation on 
the half-line. Furthermore, since f{t,0) = 99r(t, 0) = implies ip{t,0) = 0, it follows that 

/ [ipt{t,rf + (frit, rf]dr = const. (2.9) 
Jo 

The point is that the conserved quantity (j2.9p induces a stronger topology than the energy since 
it contains third derivatives of 97. Consequently, we refer to (12. 9p as a higher energy for the free 
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equation. Moreover, if we truncate the domain of integration in (I2.9p to the backward Hghtcone 
Ct, we obtain a local version of the higher energy given by 

[ [iptit,rf +ipr{t,r)^]dr. (2.10) 
Jo 

A simple scaling argument (or a straightforward computation) then shows that the higher energy 
for the blowup solution tp'^ behaves like (T — t)~^ and thus, unlike the original energy, the local 
higher energy (12.10p is strong enough to detect self-similar blowup. Consequently, we study the 
Cauchy problem Eq. (j2.7p in the topology induced by (j2.10p . 



2.3. First— order formulation. We intend to formulate Eq. (12. 7p as a first-order system in time. 
To this end, we introduce two auxiliary fields ifi, ip2 by 

Mt,r) := ^j^ipt{t,r) (2.11) 
^2{t,r) := {T-t)ldr[ldrir^^{t,r))]. 

The definition of the field ip2 is motivated by the discussion in Section [2. 2[ In fact, apart from the 
factor T — t m front, (p2 is exactly the function from Section [2.2i The factor T — t is introduced 
to put ip2 on the same scaling level as the original field ip. The field ipi is a suitably scaled time 
derivative of (p which leads to a simple expression for the higher energy (j2.10p in terms of ipi and 
ip2- Note further that, on any time slice t = const, ip can be reconstructed from ip2 by 

V^it,r) = ^^^{K\2{t,-)){r) 

where the integral operator 

Kf{p):= r sf{s)ds 
Jo 

will appear frequently in the sequel. Furthermore, a straightforward computation shows 
where K(p2 is an abbreviation for {K{p2{t, ■)){r). Consequently, Eq. (j2.7p transforms into 

2yi _L r-Va _ 3K^2 _ 3F'(V>T)-6 j^2,„ _ r ( K^'P2 

(2.12) 



9i(^2 = (T-t)3i5,(ia,(^i) 



for (fj : Ct — > ]R, i = 1, 2, with initial data 



<^i(0,r) = ^b(r)-VnO,r)] 



(2.13) 



r) = TC-drf[r\f{r) - ^"^(0, r))]. 
In order to write this in a more concise form we introduce a differential operator T)^ given by 

V^f{r) := rf"{r) + 5/'(r) + f.f{r). (2.14) 
Then we have V"^ f{r) = (-5r)^[r'^/(r)] and thus, the initial data in Eq. ()2.13p can be written as 



¥'i(0,r) = 4[5(r)-VnO,r)] 



(2.15) 



2.4. Similarity coordinates. Note that, via , both the "potential term" and the nonlinearity 
in Eq. (|2.12p depend explicitly on t. More precisely, they depend on the ratio In view of the 
self-similar character of the problem it is thus natural to introduce adapted coordinates ( "similarity 
variables") by setting 

r := -log(T-t), P:=^- 

The inverse map is given by 

i = r-e-^, r = e-^/j 

and the derivatives transform according to 

dt = e^{dr + pdp), dr = e'^dp. 

Furthermore, under the transformation {t,r) i— > {t, p), the backward lightcone Ct is mapped to the 
infinite cylinder 

Zt:={{t,p) :T>-logT, pG [0,1]}. 
Consequently, in the new coordinates (r, p) the blowup takes place at infinity. By setting 

4>j{T,p) ■■= (pj{T - e~'',e~^p), j = 1,2 
we obtain from Eq. (|2.12p the system 




pdp4>i + 201 + p^(t>2 - 3K(t>2 - V{p)K^(t>2 - pNt [j,K^2 
j = 1,2, with data 



br{-\ogT,p) = Tp%{Tp) - i^T{0,Tp)] 
b2{- logT, p) = TV^if - ^^(0, -Wp) 



(2.16) 



(2.17) 



and the potential 



Note carefully that, by transforming to similarity variables, all the explicit dependencies on the 
time variable have disappeared and we have effectively reduced the study of the self-similar blowup 
solution Tp'^ to a small data asymptotic stability problem given by Eq. ()2.16p . The analysis of 
Eq. (j2.16p is the content of the present paper. 



3. Linear perturbation theory 



In this section we study the linearized problem that results from Eq. (|2.16p by dropping the 
nonlinear term. Actually, we start with the free problem which follows from Eq. (|2.16p by dropping 
the nonlinearity and the potential term. Our approach is operator-theoretic. The point is that we 
need to employ semigroup theory in order to solve the free problem since the transformation to the 
nonorthogonal coordinate system (r, p) has in fact destroyed the underlying self-adjoint structure 
of the wave operator. Consequently, we rewrite Eq. (|2.16p as an ordinary differential equation (in 
r) on a suitable Hilbert space which is dictated by the local higher energy defined in (j2.10p . Then 
we prove well-posedness of the free problem by an application of the Lumer-Phillips theorem. The 
analogous result for the linearized problem follows by a general abstract perturbation argument, 
although the corresponding growth bound of the evolution that is obtained by this procedure is far 
from being optimal. In order to improve this bound, we perform a more detailed spectral analysis. It 
turns out that the linearized time evolution exhibits an inherent instability which is a manifestation 
of the time translation invariance of the original problem. We show how to construct a suitable 
spectral projection that removes this "artificial" instability and proceed by proving a decay bound 
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for the linearized evolution on the stable subspace. This result, which is almost optimal, concludes 
the study of the linearized problem. 

3.1. Function spaces and well posedness of the linear problem. We are going to need the 
following version of Hardy's inequality. 

Lemma 3.1. Let a > 1 and assume that u € C[0, 1] has a weak derivative as well as 

lim ' = 0. 

ThcTi 

''Hp)W^^ ( 2 V f'\u'{p)\'' 



I 

Jo 



dp < { Z 7 1 / ,a-2 '^P- 



10 Va-iy Jo P° 

Proof. This follows by integration by parts and the Cauchy-Schwarz inequality. □ 
We set 

7{ := |u= {ui,U2) e C^[0,1] X C^[0,1] : u[^\o) = 7/2(0) = 0, A; = 0,1,2,3} 
and define a sesquilinear form (•!•) by 



u V 



1 rl 

2„ 



{ui\vi)i + {u2\v2)2 ■= D ui{p)D'^vi{p)dp+ U2(p)v'2{p)dp 



where Df{p) := ^f'{p). Note that (-j-) is chosen in such a way that it leads to the local higher 
energy Eq. (j2.10p . 

Lemma 3.2. The sesquilinear form (•!•) defines an inner product on 7i and the completion ofH, 
denoted by %, is a Hilbert space. Furthermore, the subspace C^(0, 1] x C^(0, 1] ofJi is dense and 
u £ T-L implies u G C^[0, 1] x C[0, 1] with the boundary conditions tii(O) = u'^(O) = U2{0) = 0. 

Proof. From Hardy's inequality we obtain the estimate 

'\D\,{p)\'dp< f\uf{p)\^dp 
Jo 

for \i £ T-i which shows that (-j-) is well-defined on all oi H x 7i. Furthermore, the assumed 
boundary conditions ensure that (u|u) = if and only if u = 0. By the density of C^(0, 1] in 
L^(0, 1) we can, for any e > 0, find a function v € C^(0, 1] such that — 7'||l2(q ;^) < e. By 

setting V := K'^v we obtain v € C^(0, 1] with \\D'^{ui — 7^)||l2(o,i) < e which yields the claimed 
density property. Finally, we note that 

|inUp)| < \d,[^^u[{p)]\dp < \D'u,{p)\'dp 

by Cauchy-Schwarz. □ 
Now we set 

V(Lo) := |u = (Mi,n2) € C°°[0,1] x C°°[0,1] : ^(0) = ^2(0) = 0, A: = 0, 1, 2, 3, 4} 

and define a differential operator on X'(Lo) by 

-pu[{p) + 2ui{p) + p'^U2{p) - 2>Ku2{p) 

d^Mp) - pu'2{p) -Mp) 

where, as before, Kf{p) := sf{s)ds. At this point it is important to note that 

p\2{p) - SKu2{p) = 0(/) 



Lou(p) 



instead of only 0{p^) as one might expect at first glance. This is due to a special cancellation. As 
a consequence we see that 

[Unf\p) = 0{p^-^), A; = 0,1, 2, 3, 4 
where [LquJj, j = 1, 2, denotes the j-th component of Lqu. Similarly, we have 

[Lou]f (p) = 0(pi-'=) 

for A; = 0, 1 and we conclude that Lq has range in %. Comparison with Eq. (j2.16|) shows that 
Lq represents the right-hand side of the free problem and Lemma 13.21 implies that Lq is densely 
defined. Furthermore, in view of the definitions of , ip2 in (j2.1ip , the boundary conditions required 
in I'(Lo) are natural. 

Lemma 3.3. The operator Lq : V^ho) C Ti —?■ H is closable and its closure Lq generates a strongly 
continuous one-parameter semigroup Sq : [0, oo) — )• B{'H) that satisfies 

l|So(r)|| <e-^ 
for all T >0. In particular, the Cauchy problem 

r f <l>(r) = Lo<I>(t) 

\ ^'(0) = uen 

has a unique mild solution <I> : [0, oo) — )• H given by $(t) = So(t)u. 

Proof. According to the Lumer-Phillips Theorem (see |19l , p. 83, Theorem 3.15), it suffices to show 
that 

• Re(Lou|u) < -|||u|p for all u G I^(Lo) and 

• the range of A — Lq is dense in Ti for some A > — |. 

In the following, we employ a common abuse of notation and use the symbol p to denote both the 
independent variable and the identity function. Furthermore, all integrals run from to 1 and we 
omit denoting the measure dp. In order to estimate Re(Lou|u), we start by collecting all terms 
that only contain ui and integrate by parts to obtain 

-Re J D'^{p'^Dui)D^ui + 2 J jD^uip = -Re J p{D'^ui)' D'^ui - 2 j [D^uip 

= -i|I)V(l)p-i / II^VP 

where we have used the commutator [D^, p^D] = AD^ and pu'i{p) = p'^Dui{p). The boundary term 

at vanishes thanks to u^^\q) = for k = 0, 1, 2, 3, 4. Similarly, the terms containing only U2 are 
given by 

-Re j {pu2)%- j = -Re j pu'fi^ - 2 j 
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As a consequence, it suffices to show that the mixed terms are dominated by ^(|D^ui(l)p + |n2(l)P) 
and indeed we have 

Re j [D'^{p^U2) - ?,Du2]D^ui + Re j (D^ui)'^ 

= Re y M2^2^i + Re [D'^ui{1)v^T)] j {D^ui)t4 
= Re [i?V(l)^(T)] < U\D\i{l)\' + 14(1)1') 

since ^2^1(0) = 0. Thus, we obtain Re(Lou|u) < — |||u|p as desired. Note that this result is not 
surprising since it is just a reflection of the fact that the higher energy (I2.9P is conserved for the 
free problem and the factor — | can be concluded by a scaling argument. 

It remains to show that the range of A — Lq is dense in T-L for some A > — | . To this end it suffices 
to show that the equation (2 — Lo)u = f has a solution u € I^(Lo) for any f G C^(0, 1] x C^(0, 1] 
(cf. Lemma [3. 2p . The point is that the equation (2 — Lo)u = f can be solved explicitly by elementary 
ODE methods. We just state the result. For given f = (/i,/2) G C';?°(0, 1] x C^{0,1] define an 
auxiliary function u by 

n{p) := j^^^ 1^ ^ [/i(.) + s'Khis)] ds. 

Observe that u G C°°[0, 1] by Taylor expansion and u^''\p) = 0{p^-^) for k = 0,1,2,3. Now 
set U2 '■= Du. Then we have U2 G C°°[0, 1] and U2{p) = 0{p). Furthermore, define ui := 
K{p^Du) + Ku- K^f2 which implies m G C°°[0, 1] and (p) = 0{p^-^) for A: = 0, 1, . . . , 5. 
Consequently, we obtain u = {ui,U2) G ^^(Lo) and by straightforward differentiation one verifies 
that indeed (2 — Lo)u = f. Since f was arbitrary we are done. □ 

Next, we add the potential term from Eq. (j2.16p which is represented by the operator L', defined 

by 

L u(/)) := ( Q 

with the smooth potential V given explicitly in Eq. (|2.18|) . Since V G C°°[0, 1], it follows that 
L' G B{T-L) (use Hardy's inequality) and we can immediately conclude the well-posedness of the 
linearized problem. 

Corollary 3.4. The operator L := Lq + L' generates a strongly continuous one-parameter semi- 
group S : [0, 00) % which satisfies 

||S(t)|| < e(-t+ll^'ll)^ 
for all T >0. In particular, the Cauchy problem 

1^{t) = Lci>(r) 

$(0) = u£n 

has a unique mild solution <I> : [0, 00) — > 7i given by ^(t) = S(t)u. 

Proof. This is a consequence of the Bounded Perturbation Theorem, see [19j, p. 158. □ 
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3.2. Spectral analysis of the generator. In order to improve the rough growth bound for the 
Hnearized evolution given in Corollary 13.41 we have to analyze the spectrum of L. Note first that 
the growth bound for the free evolution in Lemma 13.31 implies 

cj(Lo) C {A e C : ReA < -|}, (3.19) 

see [19], p. 55, Theorem 1.10. In fact, it is not very hard to see that we have equality here, i.e., the 
growth bound in Lemma 13.31 is sharp. However, we will not need this result in the following and 
therefore we omit its proof. Of course, the addition of the potential term L' changes the spectrum; 
but as a consequence of the following result, the change is in some sense the mildest possible: it 
only affects the point spectrum. 

Lemma 3.5. The operator Jj' : Ti ^ H is compact. As a consequence, a(L)\a(LQ) C fTp(L). 

Proof. We write V. = 7ii x T-L2 and denote by |[ • j = 1,2, the respective norms on T-Lj. Since 
multiplication by V is bounded as an operator from Tii to Hi (Hardy's inequality), it suffices to 
show that is compact as an operator from 'H2 to Tii. Let (uj) C 7^2 be a bounded sequence. 
By definition of || • II2 and the boundary condition Uj{0) = it follows that (uj) C H^{0,1) is 
bounded and the compact embedding H^{0,1) CC L^(0, 1) implies that (uj) has a subsequence 
which converges in L^(0, 1). Since HiT^Ujlli = H^j ||l2(o,i) conclude that (K'^Uj) has a convergent 
subsequence in Tii which implies the compactness of K^. 

If A G (t(L)\(t(Lo) then it follows from the identity A — L = [1 — L'Rlo(A)](A — Lq) and the 
spectral theorem for compact operators (Riesz-Schauder theory, see e.g. ^T], Section 5.4) that 
A G ap(L). □ 

As we will show now. Lemma 13.51 provides the link between the spectral problem for L and the 
mode stability ODE (fLljl . 

Lemma 3.6. //AG c(L) and ReA > — | then there exists a nontrivial u G C°°[0, 1] such that 

- (1 - p") \u" + ^vl] + 2\pu' + A(A + l)u + ^^W^W) ^ ^ o_ (3 20) 

Proof. Let A G o"(L) with ReA > — |. According to Lemma 13.51 and Eq. p.l9p we have A G o'p(L) 
and thus, there exists a nontrivial u G /^(Lq) C H such that (A — L)u = 0. Writing out the 
components we obtain the two equations 

Ani(p) + pu[{p) - 2ui{p) - p^U2{p) + 3Ku2{p) + V{p)K^U2{p) = 
Au2(p) - D^ni(p) + pu'^ip) + U2{p) = 0. 

The second equation implies 

Mp) = r s^U2{s)ds + {\-l)K'^U2{p) (3.22) 
Jo 

which in particular shows that U2 is nonzero. In view of Eq. (j2.1ip we set 

u{p) := ^K^U2{p) (3.23) 
and note that u G C^[0, 1]. With this definition the expression for ui in Eq. ()3.22p simplifies to 

ui{p) = p\pu'{p) + \u{p)]. (3.24) 
Inserting Eqs. (j3.24p . (j3.23p into the first equation of (j3.2ip we infer 

-(1 - p^) \u"{p) + lu'{p)\ + 2\pu'{p) + A(A + l)u{p) + Vip)] u{p) = 



(3.21) 



p 

and by recalling the definition of V in Eq. (|2.18p we see that u indeed satisfies Eq. (j3.20p . Note 
that the coefficients in Eq. (I3.20p belong to C°°(0, 1) and furthermore, the coefficient of u" does not 
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vanish in (0, 1) which shows that the solution u is in C°°(0, 1) by basic ODE theory. The behavior 
at the endpoints follows by Frobenius' method: at /3 = the Frobenius indices are {—3,2} and 
therefore, u G C^[0, 1] already implies u G C°°[0, 1) with u{p) = 0{p'^) as p ^ 0+. At p = 1 we 
have the indices {0, 1 — A}. Note that w ^% implies U2 G H^{0, 1) and thus, u G H^{^, !)• Since 
Re(l — A) < I by assumption, the condition u G H^{^,1) excludes the nonsmooth solution at 
p = 1 and we obtain u G C°°[0, 1] as claimed. □ 

3.3. Construction of the spectral projection. As already mentioned in the introduction, the 
function g{p) := pWq{p) solves the mode stability ODE ()3.20p with A = 1. Via the transformations 
in the proof of Lemma [HTU] (in particular Eqs. (j3.22|) and (|3.23p ). g gives rise to a function g G 2?(Lo) 
which (after a convenient normalization) reads explicitly 

sip) = -m {p'[P9'{p)+gip)], j;dp[j;dp{p'9ip))]) (3.25) 

_ ( p^jb-p"^) p(125 - 50/9^ -3p4) \ 
"V(5 + 3p2)3' (5 + 3p2)4 ) 

and satisfies (1 — L)g = 0. Consequently, 1 G o-p(L) but, as already indicated in the introduction, 
this instability is induced by the time translation symmetry of the Yang-Mills equation (jl.ip . 
Our aim is to construct a suitable spectral projection that removes this symmetry mode. As a 
preparation for this we need the following observation. 

Lemma 3.7. The eigenvalue 1 G o"p(L) is isolated in the spectrum ofL and its algebraic multiplicity 
is finite. 

Proof. According to Lemma 13.61 each A G o"(L) with ReA > — | gives rise to a nontrivial function 
u G C°°[0, 1] that satisfies Eq. (|3.2Up . In fact, inspection of the proof of Lemma 13.61 shows that u is 
even analytic. Consequently, all A G cr(L) with ReA > — f are zeros of an analytic function (namely 
the Wronskian of the two analytic solutions of Eq. (j3.20p around p = and p = 1, respectively) 
and therefore they are isolated. If the algebraic multiplicity of 1 G o'p(L) were infinite then, by |21j . 
p. 239, Theorem 5.28, 1 would belong to the essential spectrum0 of L. However, since the essential 
spectrum is stable under compact perturbations ([2lj, p. 244, Theorem 5.35) and 1 ^ c7(Lo), we 
conclude that the algebraic multiplicity must be finite. □ 

Lemma [5771 allows us to define the Riesz projection 

^■■=^rli^-^r'd\ (3.26) 

where F is a circle that lies entirely in p(L) and encloses the eigenvalue 1 in such a way that no 
other spectral points of L lie inside F. By definition, the algebraic multiplicity of 1 G <7p(L) equals 
dimrgP and thus, by Lemma 13.7^ P is of finite rank. Moreover, P commutes with L in the sense 
that PL C LP and as a consequence, P also commutes with the semigroup generated by L, i.e., 
PS(t) = S(t)P for any r > 0. We set M := rgP which is a finite-dimensional subspace of Ti 
and denote by L_vi := ^'\v(L)nM the part of L in Ai. L^ is a linear bounded operator on the 
finite-dimensional Hilbert space Ai with a(Lj^) = {1}. We refer to j21j for these standard facts. 

Lemma 3.8. The subspace Ai = igP is one-dimensional and spanned by the symmetry mode g. 

^Strictly speaking, the cases A £ { — 1,0,1} require special attention since for these values of A there exist two 
possibilities: the nonsmooth solution involves a logarithmic term or all solutions are smooth at p = 1. In either case, 
however, we arrive at the same conclusion as for A ^ { — 1, 0, 1}. 

There exist at least five nonequivalent notions of essential spectra for nonself-adjoint operators, see for a 
detailed discussion. We stick to the definition given by Kato [21] as the set of all A such that A — L fails to be 
semi-Fredholm. 
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Proof. Note first that it follows from the proof of Lemma 13.61 that the geometric eigenspace of 
1 G crp(L) is one-dimensional and spanned by g. Consequently, since 1 S cr(L_A^) = crp(L_A/(), we 
conclude that g E which shows (g) C Ai. 

In order to prove the reverse implication observe that a{l — Ijm) = {0} thus, 1 — L;k is 
nilpotent. This means that there exists an m G N such that (1 — L_v()™u = for all u € 7W and 
we assume that m is minimal with this property. If m = 1 it follows that A4 C ker(l — Ijm) = (g) 
and we are done. Thus, assume m > 2. Then there exists a nonzero v € rg(l — L;k) such that 
(1— L_^)v = 0. In other words, g € rg (1— L_a/(), i.e., there exists a u € V^L) such that (1 — L)u = g. 
By a similar computation as in the proof of Lemma 13.61 we infer the equation 

-(1 - p2) + j.u'ip)] + 2pu'{p) + 2u{p) + ^I^^jPtz^uip) (3.27) 

= ^9i{p) + ^Kg,{p) - ^K^g,{p) = ^^^^ =: dip) 

for the function u{p) := ■\K'^U2{p). The homogeneous version of Eq. (|3.27p has the fundamental 
system {/io,/ii} where 

ho{p) = j^^^ 

and hi (which can also be given in closed form) behaves like \hi{p) \ ~ ^ as p ^ 0+ and |/ii(/o)| ~ 
I log(l — p)\ as p — > 1 — . After a suitable normalization of hi we obtain for the Wronskian of /iq and 
hi the expression 

W{hM{p) = ^;TlT^) 
and thus, according to the variation of constants formula, u must be of the form 

rp rp 
u{p) = coho{p) + cihi{p) + ho{p) s^hi{s)g{s)ds - hi{p) i s^ho{s)g{s)ds 

J po J pi 

for suitable constants co,ci G C and PQ,pi € [0,1]. Since limp_j.o+ ^(p) exists, we must have 

ci = f^^ s'^hQ{s)g{s)ds and thus, 

rp rp 

u{p) = cohoip) + ho{p) s^hi{s)g{s)ds - hi{p) i s'^ho{s)g{s)ds. 

J PO Jo 

Similarly, the existence of lim.p^i-u{p) yields s'^h(j{s)g{s)ds = since hi is in near p = 1. 
This, however, is impossible since s'^ho{s)g{s) > for s G (0,1). Consequently, there cannot exist 
a u € X^(L) such that (1 — L)u = g and we must have m = 1. □ 

3.4. The linear time evolution restricted to the stable subspace. As already mentioned 
several times, the unstable eigenvalue 1 G o'p(L) stems from a symmetry mode and does not 
correspond to a "real" instability. Consequently, we consider the linear time evolution on the stable 
subspace A/" := rg (1 — P) = ker P where P is the spectral projection defined in Eq. (j3.26p . Our aim 
is to derive a decay estimate for the subspace semigroup S(t)|_/^. To this end it is useful to recall 
the definition of the spectral bound of a closed operator. 

Definition 3.9. Let A : P(A) C X — > A be a closed operator on a Banach space X. Then the 
spectral bound s{A) is defined as 

s{A) := sup{ReA : A G cr{A)}. 



As before, we denote by 'Lj\f the part of L in A/" and recall that o"(L_/^) = (t(L)\{1} (see |21j). 
Note that according to numerics [2] we have in fact s(L_/^) ~ —0.59 but we emphasize that this 
information is not needed for the linear theory we are currently developing. We need a preparing 
result which is interested in its own right as it shows that there do not exist unstable eigenvalues 
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far away from the real axis. We remark that this statement does not depend on the special form 
of the potential V. It is merely a consequence of the structure of the differential operator L. 

Lemma 3.10. Set Ha := {z G C : Re z > a}, a G M. For any e > there exist constants Ci, C2 > 
such that 

I|Rl(A)|I <C2 

for all A G H_'ij^^ with |A| > Ci. In particular, L does not have unstable eigenvalues far away from 
the real axis. 

Proof. Let A G H_3_^_^ for a fixed but arbitrary e > 0. The identity A — L = [1 — L'Rlo(A)](A — Lq) 
shows that A — L is invertible if and only if 1 — L'Rlo(A) is invertible. Thus, we have to estimate 



L'R..,(A)f=(-^^'[^-"(^)f]^) 



for i £ T-L. We write u = RLo(A)f and thus, (A — Lo)u = f. The second component of this equation 
implies 



r s\2{s)ds + (A - l)K\2{p) - K^f2{p) 
Jo 



or, in other words, 

[RLo(A)f]i(/>) = r s^[KuAm2{s)ds + (A - l)K'[RuAm2{p) - K^f2{p). 
Jo 

From this we obtain the estimate 

by noting that ||Rlq(A)|| < j^_^3 (Lemma 13.31 and [19] . p. 55, Theorem 1.10) where || • j = 1, 2, 

denotes the norm on "Hj. As a consequence, if |A| is sufficiently large, the Neumann series 

00 

[1-L'Rlo(A)]-i = ^[L'Rlo(A)]'= 

k=0 

converges in norm and the claim follows. □ 

To conclude the linear perturbation theory, we estimate the linear evolution on the stable sub- 
space depending on the spectral bound of its generator. 

Proposition 3.11. Let e > and set uj := max{— |, s(L;v')} + £• Then there exists a constant 
Cg > such that the semigroup S(r) given in Corollary \3.4\ satisfies the estimate 

||S(r)(l-P)|| <C,e-^ 

for all T >0 where P is the spectral projection defined in Eq. (|3.26p . 

Proof. The operator L_/v' is the generator of the subspace semigroup S(r)|;\^ = S(r)(l — P) and its 
resolvent is given by Rl(A)|^. Consequently, the claimed estimate for S(r)(l — P) follows from the 
uniform boundedness of Rl(A) in the half-space (Lemma lS.lOp and the Gearhart-Priiss-Greiner 
Theorem ([19], p. 302, Theorem 1.11). □ 

We remark that if the self-similar solution ip'^ is mode stable (cf. Definition II. 2p then Lemma 
13.61 implies s(Lx) < 0- Hence, Proposition 13.111 shows that mode stability of ip'^ implies linear 
stability. The numerically obtained value s(L^) Rr: —0.5889 [2j yields the exponential decay 

||S(r)(l-P)|| <e-0-58.^ ^>0 

for the linearized time evolution of perturbations of ip'^ . 
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4. Nonlinear perturbation theory 



Based on Section[3]we are now ready to treat the full system Eq. (|2.16p . From now on we assume 
that s(L_/v') < 0, i.e., that ijP" is mode stable. Proposition 13. Ill then shows that the linearized time 
evolution on the stable subspace decays exponentially. This puts us in an extremely convenient 
position since normally, at least in the study of wave equations, one can at most hope for polynomial 
decay due to the continuous spectrum of the Laplacian. In fact, as is well-known from dynamical 
systems theory, exponential decay of the linearization carries over to the nonlinear evolution via 
Duhamel's formula. In the PDE context one is of course faced with the additional complication 
that one needs good mapping properties of the nonlinearity with respect to the spaces defined by 
the linear problem. However, the norm we are using controls three derivatives and we are dealing 
with a 5-dimensional problem where 1+ derivatives already suffice for a Moser estimate. It is thus 
not surprising that we are able to obtain a Lipschitz property of the nonlinearity which is necessary 
to run a fixed point argument. However, the presence of the symmetry mode g renders the linear 
evolution unstable and we have to overcome this by restricting ourselves to special initial data that 
live on a codimension one "manifold". In a second step we then remove this restriction by adjusting 
the blowup time T. At this point the role of the symmetry mode and its connection to the time 
translation invariance of the problem become evident. As a matter of fact, the symmetry mode g 
corresponds to the derivative at T = 1 of the curve T i— )• (-0^(0, ■),Tpf{0, •)) in the space of initial 
data. 

4.1. Estimates for the nonlinearity. As in the proof of Lemma 13.51 we write Ti = Hi x "^2 and 

denote by || • ||j the respective norm on Hj, j = 1,2. As a reminder we recall that 

Iklli = / \D'^u{p)\'^dp, \\u\\l = / \u'{p)\'^dp 
Jo Jo 

where Du{p) = ^u'{p). Furthermore, from now on we assume all functions to be real-valued. 

The nonlinear term in Eq. (|2.16p reads pNx{-^K'^4>2)- In fact, since we switched to similarity 

coordinates, we have by Eq. (j2.8p . 

Nt{u){p) = mo{p)u{pf + Mpf 

and thus, the nonlinear term is independent of T. For the following it is useful to separate the 
functional dependence more clearly and we therefore define 

N{x, p) := 9Wo{p)x'^ + 3x^ (4.28) 

Since the nonlinear term occurs in the first component of Eq. (j2.16p and takes an argument from 
the second component, we have to study its mapping properties as a map from to Hi. We start 
by defining an auxiliary operator 

Au(p) := ^K'^u(p) 
p-^ 

which represents the argument of the nonlinearity where, as always, Ku{p) = su{s)ds. 
Lemma 4.1. We have the bounds 

II(-)~^^«IIl2{o,i) ^ IWh 

3 

||(-)"5Au||ioo(o,i) < \\u\\2 

for all u € T-L2- 

Proof. By Lemma 13.21 we may assume u € C^(0, 1] and thus, 

ll(-) ^^lli2(oi)= / -r^dp= / J w^dp< / \u'{p)\^dp 

Jo P Jo P Jo 

16 



by repeated application of Hardy's inequality (Lemma 13. ip . For the second bound we note that 
\Au(p)\ < ^ f s f t\u(t)\dtds 

r Jo Jo 



Jo \Jo J \Jo 



^ P^||u'||i2(o,i) 



again by Hardy's inequality. □ 

We provide similar bounds for the derivatives of A. 
Lemma 4.2. We have the bounds 

II-D^^*IIl2{o,i) ^ IMI2 

||(-)'Z?'^n||,.2(o,i)<||n||2 
||(-)^L'^M||ioo(o,i) < ||n||2 

Il°°(o,i) 

for all u € 7^2- 



JL>2An||ioo(oi) < |ln||2 



Proof. The proof consists of straightforward applications of Hardy's inequality, the logic being, of 
course, that each application of D loses two powers of p. □ 

Now we define N{u){p) := pN(Au{p), p) which corresponds to the nonlinearity in Eq. ()2.16p . We 
have the following crucial result which is key to control the nonlinearity. 

Lemma 4.3. The function N maps 7i2 to Hi and we have the bound 

\\N{u)-N{v)\\i<{\\u\\2 + \\v\\2)\\u-v\\2 

for all n, ?; G B2, the open unit ball in 'H2- Furthermore, N{0) = and N is Frechet differentiate 
at with derivative D^N{0) = 0. 

Proof. Evidently, we have 

N{x, p) - N{y, p) = [QWq{p){x + y) + ?,{x^ + xy + y^)]{x - y) 

for all x, y G M and p G [0,1] and thus, we obtain 

D\N{Au) - N{Av)] = D'^[^{-)Wq{Au + Av) + ^{■){{Auf + AuAv + {Avf)]{Au - Av) 

+ 2D[Q{-)Wq{Au + Av) + ?>{-){{Auf + AuAv + {Avf)]D{Au - Av) 

+ [9{-)Wq{Au + Av) + 3(-)((An)2 + AuAv + {Av)'^)]D'^ {Au - Av). 

We need to put this whole expression in and therefore, we place the terms involving Au — Av 
in and the rest in L°° and apply Lemmas 14.11 and 14.21 to bound them. For instance, we use 

\\D\-)Au]{Au-Av)h2^o,i) < \\i-)^D\.)Au]h^^o,i)\\irHAu-Av)h2^o,i) 



< \\u\\2\\u - v\\2 



Since 

11^ \2 n2 



D2[(.)^u]||i-(0,i) < ll(-) ^^tx||i^(o,i) + ||(-)I?An||i-(o,i) 

) 



+ ||(-)'Z?'An||ioo(o,i) 



1 1 M 
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Similarly, we obtain 

< lltijll \\u — v\\2 



since, e.g. 

\\{-fD\Auf\\Lo.^o,i) < ||(-)'(^^^)'||l-(o,i) + \\{-fAuD^Au\\L^^o^^^ 

< Ml 

etc. The other terms can be treated in the exact same fashion. Since N{0,p) = for all p G [0, 1], 
we obtain A^(0) = and thus, ||A^(n)||i < ||n||2 for all u € 82- This estimate also implies the 
Prechet differentiability of TV at with D^N{0) = 0. □ 

In order to write the main equation ()2.16p as an ordinary differential equation on the Hilbert 
space 7i we introduce the vector-valued nonlinearity 

N(u) := ( -«(«^) 
for u = {ui,U2). 

Lemma 4.4. The nonlinearity N maps % to itself and satisfies the hound 

||N(u)-N(v)|| <(||u|| + ||v||)||u-v|| 

for all u, V G ;B where B denotes the open unit hall in %. Furthermore, N(0) = and N is Frechet 
differentiable at with L»^N(0) = 0. 

Proof. Since ||N(u)|| = ||A^(ti2)||i, the statement is an immediate consequence of Lemma 14. 3[ □ 

Consequently, Eq. (|2.16p can be written as 

^$(t) = L^>(r) + N(«>(r)) (4.29) 

and our aim is to study Eq. (j4.29p with small Cauchy data u G T-L prescribed at r = — logT. Thus, 
by Duhamel's formula we may rewrite the problem as 

$(t) = S(T + logr)u+ /" S(t - r')N($(r'))(ir' (4.30) 

J-logT 

which is equivalent to 

^{t) = S{t)u+ [ S{t -T')-N{^{T'))dT' (4.31) 
Jo 

for ^'(t) = $(T-logr). 

4.2. Existence for codimension one data via the Lyapunov— Perron method. Our goal 
is to prove global existence for Eq. (iOT]l . This is not straightforward since the linear time 
evolution is unstable due to the presence of the symmetry mode g given in Eq. (|3.25p . Indeed, we 
have S(T)g = e^g. Consequently, in a first step we modify Eq. (j4.3ip and consider 



roo rr 

^r(r) = S(t)(1 - P)u - / e^-^'PN {^{T'))dT' + / S(r - t')N{^{t'))t' 

Jo Jo 



(4.32) 



^Recall that global existence in the variable r really means local existence for the original equation in the backward 
lightcone Ct- 
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instead. Comparison with Eq. (j4.3ip shows that we have actuahy modified the initial data by 
subtracting the term 



/>oo 

u+ / e-^'N{^{T'))dT' 
Jo 



(4.33) 



which is an element of the unstable subspace (g). However, note carefully that the modification 
depends on the solution itself. As we will see, this modification stabilizes the evolution and we are 
able to obtain global existence. The procedure of modifying the data in order to force stability of 
the evolution is known as the Lyapunov-Perron method in (finite dimensional) dynamical systems 
theory. In a second step we then show how to obtain a solution of Eq. (j4.3ip . 
In order to be able to apply a fixed point argument, we define an operator K by 

K(^';u)(r) := S(r)(l -P)u- / e^-^'PN(^'(r'))dr' + / S(r - r')N('f (r'))r'. (4.34) 

Jo Jo 

As a consequence, fixed points of K(-;u) correspond to solutions of Eq. ()4.32p . We run the fixed 
point argument in a Banach space X defined by 

A' := G (7([0,oo),-H) : supel'^l^||^'(T)|l < 00} 

T>0 



where co is from Proposition 13.111 i.e., after fixing a small e > 0, the linear evolution satisfies 
|[S(r)(l — P)|| < e'^'^ with u) = max{— |, s(L_a/)} + e and u < hy the assumed mode stability. We 
also write 

W^Wx := supel'^l^||^(T)|| 

T>0 

for the norm on X. Furthermore, we denote hy Xs C X the closed subset defined by 

Xs:={^ eX : W^Wx < S}. 
In other words, Xs is the closed 5-ball in X. 

Lemma 4.5. Let 5 > be sufficiently small and assume ||u|| < 5^. Then the operator K maps Xs 
to itself and is contractive, i.e., 

|lK(vI/;u) -K(<I>;u)|U < - 
for all ^',<I> € Xs. As a consequence, there exists a unique fixed point o/K(-;u) in Xg. 

Proof. Note first that K(^';u) S C([0,oo),'H) for any ^ £ X and u G 7^ by the strong continuity 
of the semigroup S. We decompose the operator K according to 

K(^'; u)(r) = PK(^; u)(t) + (1 - P)K(^; u)(t). 

By using Lemma 14.41 and Proposition 13.111 we readily estimate 



|PK(^;u)(r)|| < / e"'"'||PN(^(r'))||(ir < supe2|^l"'||^(r')f / e"-(i+2|c.|)r'^^/ 

Jt t'>0 



00 



as well as 



:i - P)K(^; u)|| < e-H-||u|| + T ||s(r - r')(l - P)N(^'(T'))||(ir' 

Jo 
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and this yields K(^';u) C Xs for all £ Xs provided j|u|| < 5^. By a completely analogous 
computation we obtain the estimates 

\\PK{^; u)(r) - PK($; u)(r)|| < de-l'^l^||^ - ^p, 

11(1 - P)K(vI/;u)(t) - (1 - P)K(<D;u)(r)|| < ^e-l^l^l^ - 

which imply the claimed contraction property provided (5 > is sufficiently small. Consequently, 
the contraction mapping principle yields the existence of a unique fixed point in Xg. □ 

We obtain a global solution of the modified problem Eq. (|4.32p with small data. 

Proposition 4.6. Let U C Ti be a sufficiently small open ball with center in Ti. Then, for 
any given u £U, there exists a unique solution *(u) E Xs of Eq. (j4.32p . Furthermore, the map 
^ :U C Ti ^ X is continuous and Frechet differentiable at 0. 

Proof. The existence of ^(u) is a consequence of Lemma 14.51 Now note that 

||*(u) -*(v)|U < ||K(*(u);u) -K(*(v);u)|U + ||K(*(v); u) - K(*(v); v)||;, 

< i||*(u) - *(v)|U + ||K(*(v);u) - K(*(v); v)||;r 

for ah u, V eU hy Lemma K5\ and thus, |[*(u) - *(v)||^ < 2||K(*(v);u) - K(*(v); v)||a'. By 
definition of K and Proposition 13.111 we have 

||K(*(v); u)(t) - K(*(v); v)(r)|| = ||S(r)(l - P)(u - v)|| 

<e-l'^ln|u-v|| 

which yields ||K(^'(v);u) — K(^(v); v)||;f < ||u — v|[ and we conclude 

||*(U) - *(v)||;t < ||U- V|| 

for all u, V gU. Hence, ^ is Lipschitz continuous. Furthermore, we claim that [L'^^'(0)u](t) = 
S(t)(1 — P)u. Indeed, since *(0) = 0, we obtain 

e^-^'PN(*(u)(T'))(ir' 

(1 - P)[*(u)(r) - *(0)(t) - S(r)(l - P)u] = T S(t - t')(1 - P)N(*(u)(r'))cZT' 

Jo 

and as in the proof of Lemma 14.51 this implies 

||*(u)(t) - *(0)(t) - S(r)(l - P)u|I < e-H-||*(u)|a < e-H-||uf 
by the above and the claim follows. □ 

4.3. Existence for general small data. Li this section we construct a global solution to Eq. (j4.3ip . 
This amounts to removing the modification which led from Eq. ()4.3ip to Eq. (|4.32p . Recall that 
we had to introduce this modification because of the instability of the linear evolution and this 
instability emerges from the time translation invariance of the original equation. Consequently, one 
should be able to remove the instability by shifting the blowup time. Mathematically, this man- 
ifests itself in the fact that the modification turns out to be identically zero once we have chosen 
the correct blowup time T. 

Since the translated equation ()4.3ip is independent of T, the only place where the blowup time 
enters is in the data. For given data (/, g) as in Eq. (12. 6p . we denote by v the corresponding data in 
the new coordinates relative to ip^, the fundamental self-similar solution with blowup time T = 1. 
Explicitly, we have 

V^^IZ-V^Ho, •)](/>) j' ^^-^^^ 
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cf. Eq. (I2TTD and recall that V^fir) = rf"{r) + 5/'(r) + f /(r). Then we set 

U(v,J . - ^ ^^^^^^^ J + ^ rp2[v,i(o,.) - V^(0,.)](Tp) 



(4.36) 

Tp^[giTp)-i^f{0,Tp)]^ 



Thus, with Eq. ()4.35p . the initial data for Eq. (|4.29p can be written as 

^(0) = c^(-logr) = U(v,T), 

see Eq. (j2.17p . The point of this notation is, of course, that v is independent of T and thus, the 
functional dependence of the data ^'(0) on {f,g) and T is now explicit. We also remark that the 
data (/, g) have to be prescribed on the interval [0, T] but T is not known in advance. However, 
this defect is easily remedied by simply prescribing the data on [0, |] since we may always assume 
that T G / := (i, I) by the perturbative character of our construction. Consequently, we set 



^2' 2 

l|u||^ 



\D\i{p)\'dp+ / \u'^{p)\'dp 



and denote the respective Hilbert space by 'H. Note carefully that in view of Eq. ()4.35p we have 

l|v||^ = ||(/,5)-(^'(0,-),V't'(0,-))ll£(|) 
with the £^(|)-norm from Theorem 1 1.3i 

Lemma 4.7. The function U maps Ti x I to H continuously and U(0, 1) = 0. Furthermore, 
U(0, ■) : I H is Frechet differentiable with partial derivative 

L>fU(0,l)A = -240Ag 

for all X where g is the symmetry mode given in Eq. (|3.25p . 

Proof. It follows immediately from the definition of U that 

||u(v,r)-u(w,r)|| <||v-w||^ 

for all v,w £ 7i, uniformly in T G /. Thus, in order to show continuity of U, it suffices to prove 
continuity of U(v, ■) : I ^71 iox fixed v € 7^. By Lemma [3^ we may assume Vj S C[0, |], j = 1, 2. 
Furthermore, recah that ■0"^(O,r) = Woly) - 1 and i/j'[{0,r) = ^Wq(^) with 

This implies 'D'^ip^ {0, ■) G C°°[0, |] and continuity of U(v, •) follows by the continuity of T i— > 
|[/(T-)||l2(o,i) : / ^ M for / G C[0, |]. Obviously, we have U(0, 1) = and the Frechet differentia- 
bility of U(0, •) on / is also evident. By straightforward differentiation we obtain 02^(0, 1)A = 
— 240Ag as claimed. □ 

By Lemma 14.71 we infer that U(v, T) H is small provided v is sufficiently small in Ti and T 
is sufficiently close to 1. Consequently, we obtain U(v,r) G W where U is from Proposition 14.61 
and there exists a solution ^(U{v,T)) G of Eq. ()4.32p with initial data U(v,T). The correction 
term (I4.33p . which was introduced to suppress the instability of the linear evolution, is given by 



U(v,r)+/ e-"'N(*(U(v,r))(r'))(ir' 
Jo 



=:F(v,r) 
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and by the above considerations, F : V x J — )• (g) is a well-defined map for V a sufficiently small 
open ball around in 7^ and J C / a sufficiently small open interval with 1 G J. If F(v, T) = then 
the correction term vanishes and ^(U(v,T)) is also a solution to the original equation Eq. (I4.3ip . 
Obviously, we have F(0, 1) = since U(0, 1) = and the corresponding solution is *(0) = by 
the uniqueness in (Proposition 14. 6| ). Now we show that for any small v we can find a T such 
that F(v,T) = 0. We need one additional technical result. 

Lemma 4.8. The mapping F : V x J ^ (g) is continuous. Furthermore, F(0, •) : J — ?> (g) is 

Frechet differentiable at 1 with derivative 

L>fF(0,l)A = -240Ag 

for all A G M. 

Proof. It is convenient to introduce a symbol for the integral operator in the definition of F and 
we write 



poo 

:= / e'^'^{T')dT'. 
Jo 



Then B : A" — )• 7^ is linear and bounded. Furthermore, we define N : <Y — )■ by N(^')(r) := 
N(^'(t)). By Lemma Ol we have 

IINWIU =supeH-||N(M/(T))|| <supeH-||.I/(r)f < 

r>0 T>0 

which shows that N is Frechet differentiable at with N(0) = and D^N(O) = 0. Consequently, 
the function F can be written as 

F(v, T) = P [U(v, T) + BN(*(U(v, T))) 

and by Lemmas 14. 4^ 14.71 and Proposition 14.61 it follows that F is continuous. Furthermore, by the 
chain rule for Frechet derivatives we immediately infer 

F(0, 1) = PDf U(0, 1) + BD^N{0)D^^{0)D^ U(0, 1) 

= PD^\J{0, 1) 

and Lemma [121 yields D^F{0, 1)A = -240Ag for ah A G M. □ 

Lemma 4.9. Let V C Ti be a sufficiently small open ball around 0. Then, for any v G V, there 
exists a T G |) such that F(v,r) = 0. 

Proof. Denote by i : (g) — ?• M the vector space isomorphism given by «(Ag) = A, A G M, and set 
/ := ioF. We have /(0, 1) = and Lemma 148) shows that /(O, •) : J — > R is differentiable at 1 with 
d2f{0, 1) 7^ 0. Consequently, we obtain T^,T^ G J such that /(O, T_) < and f{0,T^) > 0. Since 
/ : V X J — > M is continuous, we find that /(v, T_) < and /(v, T+) > for all v G V provided 
V is sufficiently small. Thus, by the intermediate value theorem there exists a T G (T_,T+) such 
that /(v,T) = 0. □ 



We formulate the main result theorem. 
Theorem 4.10. Let V dT-L be a sufficiently small ball with center 0. Then, for any v G V, there 

'I 3> 
>2' 2' 



exists a T ^ (^) i) such that the Cauchy problem 



|:$(r) = L$(r) + N($(r)) 
<I>(-logT) =U(v,T) 

has a unique mild solution <I> G C{[—\ogT,oo),'H) satisfying 

II^WII <e-H- 
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for all T > -logT. 



Proof. It only remains to prove that the solution is unique in C([— logT, oo),T-L). However, this is 
a simple consequence of the fact that $ is a fixed point of a contraction mapping. □ 
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